Relevance of the slowly-varying electron gas to atoms, molecules, and solids 
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Electron densities become both large and slowly- varying in a certain semi-classical limit: The local 
density approximation becomes exact for both the kinetic and exchange energies, while the leading 
correction is given by the second-order gradient expansion. But real matter has nuclear cusps, which 
produce "Scott corrections". Thus gradient expansions must be generalized even for neutral atoms 
of large atomic number, but these generalizations should recover the correct second-order limits. 
Correlation is more subtle. 



PACS numbers: 

In a remarkable series of papers towards the end of 
his life, Schwinger[l] (sometimes with Englcrt[2]) put the 
semiclassical theory of neutral atoms on a firm footing. 
They carefully proved a variety of results, including a 
clear-cut demonstration that the local density approxi- 
mation (LDA) becomes exact for exchange as Z, the nu- 
clear charge, tends to oo. The large- Z expansion of the 
energy of atoms whose rigor they established: 

E = -0.7687 Z 7/3 + 0.5 Z 2 - 0.2699 Z 5/3 + ... (1) 

is extremely close to the total (Hartree-Fock) energies of 
neutral atoms (less than 0.5% error for Ne), an example 
of "the principle of unreasonable utility of asymptotic 
estimates." [1] But their derivations are specific to neutral 
atoms, and they eschew exploring any relation to "the 
density-functional formalism," [2] preferring to express all 
quantities in terms of the potential. Their results are 
little used within density functional theory (DFT). 

In the quarter century since, Kohn-Sham DFT has 
become a widely-used tool for electronic structure cal- 
culations of atoms, molecules, and solids [3]. Here, the 
non-interacting kinetic energy, T s , is treated exactly, and 
only the density functional for the exchange-correlation 
energy, E xc [n], must be approximated. For E xc , LDA[4] 
is very reliable, but insufficiently accurate for chemical 
bonding. The gradient expansion, employing density 
gradients, can be derived from the slowly-varying elec- 
tron gas for small gradients, but fails for real molecules 
and solids, whose gradients are not small. The devel- 
opment of modern generalized gradient approximations 
(GGA)[5, 6], depending on density gradients beyond 
that of leading order, improved accuracy and led to the 
widespread use of DFT in many fields. Anomalously, 
modern successful GGA's for exchange have gradient co- 
efficients that are about double that of the gradient ex- 
pansion, and the relevance of even LDA to exponentially 
localized densities is often questioned [7]. 

But Thomas- Fermi (TF) theory [8, 9], with its local 
approximation to T s , is the simplest, original form of 



DFT, and yields the leading term in Eq. (1). Since 
E = —T for atoms and correlation is O(Z), Eq. (1) is 
an expansion for T s . The Scott correction [10] {Z 2 term) 
arises from the cusp at the nucleus and the ls-region 
electrons [1], while the Z 5 ^ 3 term includes second-order 
gradient contributions to T s . 

We introduce a methodology that generalizes 
Schwinger's results to all systems and to other compo- 
nents of the energy. It explains why local approximations 
become exact for large numbers of electrons, and when 
the gradient expansion is accurate for real matter. It 
explains the doubling of the coefficient for exchange, 
yielding a nonempirical derivation of the fit parameter 
of the B88 functional[5], and why the analog of Eq.(l) 
fails for correlation. 

To begin, define the scaled density 

n c (r) = C 2 n(C 1/3 r), < C < oo (2) 

for any electron density. If n(r) contains N electrons, 
n^(r) contains CN electrons. We only ever consider in- 
teger (N, but study energies as smooth functions of £. 
Changing £ is exactly equivalent to changing Z in TF 
theory, and approximately so in reality. This scaling is 
defined for all systems, not just atoms. 

In Fig. 1, we illustrate the effect of the scaling on the 
density of the He atom. The density at ( — 2 contains 4 
electrons, has the original shape, but is coordinate-scaled 
by 2 1 / 3 . As C grows, the density becomes both large 
and slowly-varying on the scale of the local Fermi wave- 
length, A F (r) = 27r/fc F (r), where k F (r) = (Stt 2 ?^!-)) 1 / 3 ). 
Under this scaling, the dimensionless gradients s = 
| Vn\/(2kFn) and q = V 2 n/(4fcf.n)[4, 11] vary as 

Sc (r) = s(C 1/3 r)/C 1/3 , q C (r) = q(( 1/3 r)/( 2 / 3 . (3) 

In Fig. 2, we plot s for accurate densities (self-consistent 
OEP densities with exact exchange) of both Kr and Rn. 
The gradient expansions for the kinetic and exchange en- 
ergies are expected [4] to become exact as the density 
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FIG. 1: Scaled radial density of He atom. 
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FIG. 2: Reduced density gradient for noble gas atoms. 

n(r) becomes slowly-varying on the scale of A^, i.e., when 
s — > and \q/s\ — > everywhere. The second condition 
is needed because s — > for an infinitesimal-amplitude 
but rapid variation around a uniform density. 

Under ^-scaling of a mono- or polyatomic density, as 
£ — > oo and as s and |?/s| — > (except in regions to be 
discussed), the gradient expansion becomes asymptoti- 
cally exact so that (aside from possible oscillations [2]): 

T s [n c ] = C 7/3 T s (0) [n] + C 5/3 Ts (2) N+CT s (4) [n] + ..., 
£ x [n c ] = C 5/3 4 0) N + C4 2) [n] + ..., (4) 

where [n] is the j-th order contribution to the gradi- 
ent expansion of T s (i.e., is the TF kinetic energy, 

(2) 

Ts is 1/9 the von Weizsacker term[12], etc.), and sim- 
ilarly for exchange[13]. The terms displayed in Eq. (4) 
are those of the gradient expansion that remain finite for 
exponentially localized densities, and also those for which 
the gradient expansion becomes asymptotically exact as 
£ — > oo for analytic densities (unless the density has nu- 
clear cusps, when must be excluded). The gradient 
expansions of Eq. (4) are "statistical" approximations in 
the sense that they become relatively exact as (N — > oo. 



Evanescent regions are classically-forbidden 
regions in which the kinetic energy density 
T '(r) = £^*(r)* (-V 2 /2) Vi(r) of the Kohn-Sham 
orbitals is negative and the gradient expansion must 
fail. To investigate the contributions to T s [n^] and 
E x [n^] from evanescent and nuclear cusp regions, we 
note that each can be represented by an exponential 
density aexp(— br). For this density, the evanescent 
region r > r e (£) is defined loosely by the conditions 
s > 1 and \q/s\ > 1 (found from the second-order 
gradient expansion of r'(r)). Its radius r e (C) decreases 
like ln(^)/^ 1 / 3 as £ — > oo. In this region, the considered 
terms of the gradient expansion each contribute of order 
£ 2 / 3 to T s , and of order C 1 / 3 to E x (with the convention 
that ln(£) is of order 1). All of these contributions are 
of lower order than those shown in Eq. (4), and thus 
asymptotically unimportant. Similarly, the nuclear cusp 
regions r < r c (£) can be defined loosely by the condition 
\q/s\ > 1; their radii r c (£) decrease like £~ 2 / 3 as £ — > oo. 
In these regions, the local contributions dominate. The 
local terms contribute of order £ 4 / 3 to T s and £ 2 / 3 to 
E x , and are again asymptotically unimportant (once we 
truncate the gradient expansions of Eq. (4)) at second 
order). Our evanescent regions coincide with the edge 
surface of Ref. [14], but the cusps do not. 

While the total electron number scales up as £, the 
number in the evanescent and cusp regions remains of 
order £°. The exact kinetic energy contributions of these 
regions are of order £ 2 / 3 and £ (or ( 2 ^ 3 , for the positive 
kinetic energy density r(r) = ^ | V^i(r) | 2 /2) , respec- 
tively. All orders of £ are unchanged if the bound 1 is 
replaced by a smaller positive value, or if the spherical 
evanescent region is replaced by a planar one. 

To appreciate the significance of asymptotic exactness, 
£-scale any density (even a single exponential). Using 
modern linear response techniques, construct the KS po- 
tential and orbitals for each value of £[15]. As £ — > oo, 
deduce the local and gradient expansion approximations 
for T s and E x exactly, without ever studying the prop- 
erties of the uniform or slowly- varying electron gas. The 
system is becoming increasingly semi-classical, and the 
orbitals well-described by their WKB approximations. 
Semi-classical approximations require only that the po- 
tential be smooth locally, not globally. 

The behavior of correlation under this scaling is how- 
ever recalcitrant. In the bulk of the density, the gradient 
relative to the screening length, t — |Vn|/(2fc s n), where 
k s (r) = y/4kp(r)/n, controls gradient corrections to cor- 
relation. But 

i c (r)=i(£V3 r ) (5 ) 

does not change under our scaling, so that the density 
does not become slowly- varying for correlation. Nonethe- 
less, the local approximation still becomes exact, if the 
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PBE GGA is a reliable guide (as motivated later): 

E c [n c ] = A c (\n( + B c ( + ... (6) 

where A c = -0.02072 is correctly given by LDA, while 
£ LDA = -0.00452 and B^ BE = 0.03936. These num- 
bers can be found by applying these approximations to 
the TF density. Because t never gets small, B c has not 
only LDA and GEA (gradient expansion to second order) 
contributions, but higher-order contributions too. This 
expansion is far more slowly converging than those of T s 
and E x , and much less relevant to real systems. 

TABLE I: Noble gas atomic XC energies compared with local 
and gradient expansion approximations (hartree). 





E x 


Ec 


atom 


LDA 


GEA 


exact 


LDA 


GEA 


exact 


He 


-0.884 


-1.007 


-1.026 


-0.113 


0.103 


-0.042 


Ne 


-11.03 


-11.77 


-12.10 


-0.746 


0.559 


-0.390 


Ar 


-27.86 


-29.29 


-30.17 


-1.431 


1.09 


-0.722 


Kr 


-88.62 


-91.65 


-93.83 


-3.284 


2.06 




Xe 


-170.6 


-175.3 


-179.1 


-5.199 


3.15 




Rn 


-373.0 


-380.8 


-387.4 


-9.026 


4.78 





In Table I, we list the exact [16], LDA, and GEA results 
for the noble gas atoms using the known second-order 
gradient terms[13, 17]. LDA gets relatively better and 
better as Z grows, and the GEA for exchange does even 
better, but the GEA for correlation strongly overcorrects 
the LDA, making correlation energies positivc[17]. Thus 
this simple analysis explains why the gradient expansion 
yields good answers for T s and E x , but bad ones for cor- 
relation. These explanations complement those already 
existing based on holes and sum- rules [18], but add the 
crucial ingredient that under this scaling the gradient 
expansion becomes asymptotically exact for T s and E x . 

But Eq. (4) for cither E x or T s is far less accurate at 
( = 1 than Eq. (1) is for E at Z = 1. Careful inspec- 
tion of the origin in Fig. 2 shows that the exact curves 
approach a finite s-value at the origin, about 0.376, the 
hydrogenic value. The large- Z expansion is not the same 
as scaling to large £, except in TF theory. The gradients 
near the nucleus do not become small on the Fermi wave- 
length scale, no matter how large Z is. This region will 
contribute a term of order Z 2 to the kinetic energy and 
of order Z to the exchange energy at all levels, from LDA 
to exact. LDA applied to the TF density produces the 
leading term in E x , 0.2208Z 5 / 3 . Gradient corrections 
are of order Z, but so too is the cusp correction, i.e., 
the asymptotic expansion in large Z inextricably mixes 
these contributions (unlike in T s ). Table I lists E x for 
noble gas atoms, calculated at the self-consistent non- 
relativistic OEP level. We fit (E x + 0.2208 Z 5 ' 3 )/Z as a 
function of Z -1 / 3 , finding: 

E X (Z) = -0.2208 Z 5/3 + (C LDA + AC) Z + ... (7) 



Extraction of C LDA from LDA energies is difficult, be- 
cause of shell-structure oscillations. We estimate > 
C LDA > —0.03. However, for any other calculation of 
E x , we find (E x - E X BA )/Z is smooth, with fit results 
shown in Table II. The large underestimate of GEA shows 
that gradient corrections from the slowly-varying gas ac- 
count for only half the entire contribution. Assuming 

C LDA = 0, E X (Z) 0.2208 Z 5 / 3 - .196 Z. This yields 

less than 10% error for He, and less than 2% for Ne. 

TABLE II: AG = lim z ^ 00 (S x - Eh^ A )/Z (hartree). The 
"gradient" contribution arises from the expansion to order 
V 2 , while the remainder is the "cusp". 



AC 


GEA PBE B88 TPSS exact 


total 


-.098 -.174 -.202 -.159 -.196 


gradient 
cusp 


-.098 -.174 -.217 -.0977 -.0977 
0.00 -.000 +.015 -.0617 -.0979 



Popular GGA's such as PBE and B88[5] have second- 
order gradient coefficients that are about twice the cor- 
rect coefficient for a slowly- varying density, as they must 
to reproduce accurate exchange energies of atoms. But 
with an incorrect coefficient, they cannot predict accu- 
rate surface energies for metals [19]. The origin of the 
enhanced gradient coefficient of the GGA for exchange 
now has a simple explanation. In order to be asymptoti- 
cally exact for large Z, and hence accurate for most finite 
Z, the functional accounts for both the slowly-varying 
term and the cusp correction. No GGA can get both 
effects right individually. B88 is closest to being exact 
for AC, because of the fitting to noble gas atoms [5]. 
(In fact, assuming the exact gradient and cusp contri- 
butions to AC are equal, asymptotic exactness requires 
(3 = 5/(108(6tt 5 ) 1 / 3 ) in B88, close to the fitted value of 
.0042.) PBE preserved the nearly correct uniform-gas 
linear response of LDA for XC together [6], which pro- 
duces a AC also close to exact, and a GGA close to that 
of a hole modcl[18]. A "buried Is" region has small s, 
looking to a GGA like a region of slowly-varying density. 

However, a meta-GGA that employs -r(r) can recog- 
nize that this is a rapidly- varying region, and thus get 
everything right. The TPSS meta-GGA recovers the gra- 
dient expansion to fourth order [20], while yielding a good 
estimate for AC (Table II). 

Figure 3 plots the difference in exchange energy den- 
sities relative to LDA in the different approximations. 
PBE simply mimics GEA, being almost a factor of 2 
larger everywhere. But TPSS produces a much greater 
contribution from the region near the nucleus (via Fig. 
1 of Ref. [20]), while reverting to the GEA value at the 
inner radii of the other atomic shells, where s is small. 

Again, correlation is less clear-cut. In Fig. 4, we plot 
the correlation energy in three different approximations: 
LDA, PBE, and Moller-Plesset second order perturbation 
theory (MP2)[21]. (TPSS is almost identical on this scale 
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FIG. 3: Radial exchange energy density differences for Xe for 
different approximations (atomic units). 
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FIG. 4: Scaled correlation energies for noble gas atoms 
(hartree).The dashed lines correspond to the high-density 
limit of Eq. (6) for LDA and PBE. 



to PBE.) Also included are dashed lines that correspond 
to the high-density limit of Eq. (6) using LDA and PBE 
inputs for B c . Real atoms are so far from the asymptotic 
limit for correlation that asymptotic exactness is much 
less relevant in this case. 



TABLE III: Noble gas atomic correlation energies (hartree). 
The fixed-node Diffusion Monte Carlo (DMC) values[22] are 
upper bounds. 



atom 


LDA 


PBE 


TPSS 


MP2 


DMC 


exact 


He 


-0.113 


-0.042 


-0.043 




-0.042 


-0.042 


Ne 


-0.743 


-0.351 


-0.354 


-0.388 


-0.376 


-0.390 


Ar 


-1.424 


-0.707 


-0.711 


-0.709 


-0.667 


-0.722 


Kr 


-3.269 


-1.767 


-1.771 


-1.890 


-1.688 




Xe 


-5.177 


-2.918 


-2.920 


-3.089 


-2.647 




Rn 


-9.026 


-5.325 


-5.33 


-5.745 







Table III shows the correlation energies of closed-shell 
atoms as predicted by PBE, TPSS, and MP2 [21], along 



with essentially exact values [16] where known. The 
agreement among these values is generally good. The 
large- Z limit is not problematic for PBE and TPSS, since 
i(r) remains bounded except near the cusp and in the tail. 

Last, we relate this scaling to others. The most stan- 
dard is uniform coordinate scaling [23] (7 s 71(71-)), under 
which T s /7 2 and -E x /7 remain unchanged. More recently, 
number-scaling, in which n(r) becomes vn{r), has been 
proposed [24]. In the large v limit, all gradients in the 
bulk become small on both local length scales, making 
even the gradient expansion for correlation asymptoti- 
cally exact. The present scaling can be regarded as a 
product of these, with v = 7 s = £. A slowly-varying 
product with bounded density is n(z/ _1 / 3 r), v — > 00. 

Our ^-scaling allows the results of Schwinger's deriva- 
tions to be applied throughout DFT, yielding insight into 
the performance of approximate functionals. Even for 
uncondensed matter, such functionals should incorporate 
the second-order gradient expansions (although GGA to- 
tal exchange energies then degrade due to Is regions). 
Based on this, we have developed such a GGA for ex- 
change which is currently being tested. We thank Eber- 
hard Engel for the use of his atomic OPMKS code, and 
NSF (CHE-0355405 and DMR-0501588) and the Norwe- 
gian Research Council (148960/432) for support. 
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